Single-ion anisotropy in the gadolinium pyrochlores studied by an electron 

paramagnetic resonance 



V. N. Glazkov, 1 ' 2 M. E. Zhitomirsky, 1 A. I. Smirnov, 2 
H.-A. Krug von Nidda, 3 A. Loidl, 3 C. Marin. 1 and J.-P. Sanchez 1 

1 Commissariat a I'Energie Atomique, DSM/DRFMC/SPSMS, 38054 Grenoble, Cedex 9, France 
2 P. L. Kapitza Institute for Physical Problems RAS, 111334 Moscow, Russia 
3 Experimentalphysik V, EKM, Institut fur Physik, 
Umversitat Augsburg, 86135 Augsburg, Germany 
(Dated: February 2, 2008) 

The electron paramagnetic resonance is used to measure the single-ion anisotropy of Gd 3+ ions in 
the pyrochlore structure of (Yi-^Gd^T^O? . A rather strong easy-plane type anisotropy is found. 
The anisotropy constant D is comparable to the exchange integral J in the prototype Gd2Ti2C>7, 
D ~ 0.75 J, and exceeds the dipolar energy scale. Physical implications of an easy-plane anisotropy 
for a pyrochlore antiferromagnet are considered. We calculate the magnetization curves at T — 
and discuss phase transitions in magnetic field. 
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Pyrochlore antiferromagnets have been actively stud- 
ied in the past due to their unusual properties i The 
magnetic ions in these compounds form a network of 
corner-sharing tctrahedra, which is prone to a high de- 
gree of geometric frustration. The ground state of a clas- 
sical Heisenberg antiferromagnet on a pyrochlore lattice 
is macroscopically degenerate and remains disordered^ 
Since weaker residual interactions are always present in 
real magnetic materials, an important question is about 
a separation of energy scales. The corresponding knowl- 
edge helps to clarify how the macroscopic entropy freezes 
as T — ► and what type of magnetically ordered, spin- 
liquid, or glassy state is eventually stabilized. 

Gadolinium pyrochlores Gd2Ti20y and Gd2Sn207 are 
considered to be the best realizations of the Heisenberg 
model on a pyrochlore lattice, since Gd 3+ is, nominally, 
in a 8 <SV/2 state with completely frozen orbital degrees 
of freedom. Both compounds order at temperatures of 
about 1 Kt&ii£ while Gd2Ti20y exhibits also multiple 
phases in magnetic fieldi^ This has been attributed 
to weaker dipolar interactions between the S = 7/2 
gadolinium spins j&L£ though no unambiguous explana- 
tion neither of the phase diagram of Gd2Ti20y nor of the 
neutron-diffraction measurements^ has been presented. 

Single-ion effects have so far been neglected for the two 
gadolinium pyrochlores. However, strong spin-orbit cou- 
pling of the 4/ electrons breaks the LS-scheme of the 
energy levels of the Gd 3+ ion and mixes the 8 £7/2 and 
6 Pf/2 states^ Since the ground-state multiplet contains 
an admixture of the 1^0 states, it can be split by the 
crystal field. Crystal-field splitting of the order of 0.8 K 
was observed for Gd 3+ ions in yttrium-gallium garnet ™i 
The amplitude of such a splitting is comparable to the ex- 
change coupling and the dipolar energy in the prototype 
Gd 3 GasOi2. Thus, an accurate study of the single-ion 
anisotropy in the Gd-based pyrochlores is necessary to 
understand the magnetic properties of these materials. 

To study properties of isolated Gd 3+ ions we have pre- 



pared single-crystals of non-magnetic Y2Ti207 diluted 
with a small amount (nominally 0.5%) of Gd. The crys- 
tals were grown by traveling solvent floating zone tech- 
nique. After growth the oriented single-crystals were an- 
nealed in an oxygen atmosphere. The absence of para- 
sitic phases has been verified by powder X-ray diffrac- 
tion. The magnetization of the reference sample of 
nominally pure Y2Ti207 corresponds to no more than 
0.06% of S = 1/2 impurities per yttrium. The lat- 
tice parameters and geometry of the oxygen surround- 
ing of the rare-earth site are very close for Y2Ti207 and 
Gd2Ti207ii£ Hence, single-ion parameters for Gd 3+ ions 
in (Yo.995Gdo.oo5)2Ti 2 07 and Gd 2 Ti 2 7 are also ex- 
pected to be very close. An electric field gradient at 
the Gd site in Gd2Sn207 determined from the quadrupo- 
lar Mossbauer spectra™ is 30% smaller than that of 
Gd2Ti207, which should lead to the corresponding de- 
crease of the single-ion anisotropy. 

Electron paramagnetic resonance (EPR) is a conven- 
tional tool to probe small crystal-field splitting of the 
energy levels of magnetic ionsii 3 - The zero-field splitting 
of the multiplet leads to an easily detected change of the 
resonant absorption frequencies. We have used a set of 
home-made spectrometers with transmission type cavi- 
ties for the high-frequency (20-115 GHz) measurements, 
and a commercial Brukcr EPR spectrometer for the high- 
precision X-band (9.36 GHz) and Q-band (34.0 GHz) 
measurements . 

The local surrounding of the rare-earth site in the py- 
rochlore structure corresponds to a trigonal point sym- 
metry D^d (3m). A general form of the single-site Hamil- 
tonian is, then, 

n = b16\ + b\&{ + b\6\ + b°6° 6 

+ BlOl + BlOl -g f i B H-S , (1) 

where the Stevens operatorsi 3 *^ Oj are functions of the 
components of the total angular momentum S = 7/2, 
with z || (111) and x || (112). The lowest-order 
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FIG. 1: Observed derivative of the 9.36 GHz microwave ab- 
sorption (solid); modeled curve best fitting to the experimen- 
tal results (dashed) ; parts of the modeled curves (dotted) cor- 
responding to the spin-Hamiltonian constricted to the second- 
order term only (left) and second- and fourth-order diagonal 
terms (right). f=9.36 GHz, H||[110], T=6.0 K. 



anisotropy term -B^O" can be rewritten up to an additive 
constant as DS 2 Z with D = 3B%. 

In the absence of the single-ion anisotropy Hamilto- 
nian contains only the Zeeman term. In this case the 
frequencies of all allowed transitions are equal and only a 
single resonance absorption line should be observed at a 
field determined by hw = g^sH. The measured EPR ab- 
sorption spectra, see Figs.nand|21 contain multiple com- 
ponents demonstrating the presence of a non-negligible 
anisotropy. The analysis of the EPR spectra is compli- 
cated by the existence of four incquivalent magnetic ion 
positions. At high frequencies and H || [111] it is possi- 
ble to separate seven absorption components, which cor- 
respond to the spins with anisotropy axis parallel to the 
magnetic field direction. They include absorption com- 
ponent with the largest value of the resonance field and 
are almost equidistant (see Fig. |2J) . This indicates that 
the second order axial term (DS 2 ) is the most important 
one, its magnitude estimated from the total splitting of 
the spectra is \D\ ~ 0.25K. In order to refine the values of 
the anisotropy constants, we have modeled numerically 
the EPR absorption^ for different directions of magnetic 
field and different microwave frequencies trying to obtain 
best correspondence of the observed and modeled reso- 
nance fields values (see, e.g., Fig. To determine the 
sign of the main anisotropy constant (resonance-field po- 
sition determine only relative signs of the anisotropy con- 
stants), we have studied the temperature dependence of 
the intensity of different components at high microwave 
frequency. This analysis yielded that the main anisotropy 
term is of the easy-plane type and the spin Hamilto- 
nian parameters are g = 1.987, B\ = (74.3 ± 0.3)10" 3 K, 
B\ = (5.7±0.9)1CT 5 K, B\ = (1.8 ± 0.2)1CT 3 K. The de- 
termination of the 6-th order anisotropy parameters has 
been impossible due to a presence of small (0.5° or less) 



FIG. 2: Left panel: Field dependence of the microwave ab- 
sorption for the high-frequency EPR spectra at T—4.2K and 
H||[lll]. Dotted lines (/ = 32.07GHz) mark positions of the 
absorption determined by the spins with the anisotropy axis 
parallel to the field direction. Right panel: Field dependence 
of the selected resonance frequencies for the spins with the 
anisotropy axis parallel to the field direction ([111]). Solid 
lines give theoretical fit. 



misorientation of the sample. We were able to estimate 
only the upper bounds: < 10" 6 K, < 2 • 10~ 5 K, 
\B%\ < 10" 6 K. 

The measured value of the B% coefficient indicates that 
the single-ion anisotropy is the second largest interac- 
tion in the effective spin Hamiltonian of Gd2Ti2C>7. In- 
deed, the nearest-neighbor exchange integral is estimated 
from the Curie-Weiss temperature 3 *^ dew = —9.5 K 
as JS 2 w 3.7 K, whereas DS 2 w 2.73 K. The dipo- 
lar constant Edd = {g^sS) 2 j (a^/2/4:) 3 for the interac- 
tion between adjacent magnetic ions for above values 
of the g-factor and lattice parameter amounts to only 
Edd ~ 0.65 K. In the following, we discuss several prop- 
erties of magnetic pyrochlore materials, which might be 
affected by a strong single-ion anisotropy of the easy- 
plane type. 

The first question to be addressed theoretically is 
whether a significant single-ion anisotropy modifies the 
value of the exchange constant obtained from 0cw- Stan- 
dard high-temperature expansion^ of the susceptibility 
tensor includes in general a contribution from the single- 
ion term DJ2ii n i ' Si) 2 : 



a/3 



g 2 y 2 B S{S+l) 



3T 
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3T 



-f£(s-£)(s + §)(<^>-^)}. (2) 

In the case of a pyrochlore lattice the averaging in Eq. (J2J) 
has to be done over four inequivalent directions of the 
local anisotropy axis corresponding to the principal 

cube diagonals. Using the identity \ Ylt=i n ? n i = l^™' 3 
we find that the single-ion contribution vanishes in the 
second-order of the 1/T expansion. Therefore, the Curie- 
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Weiss constant is not affected by the single-ion term and 
is given by the same expression as in a pure Heisenberg 
case: 6cw = zJS(S + l)/3, where z = 6 is the number 
of nearest-neighbor magnetic atoms. 

The unusual ordering wave- vector (555) in Gd2Ti2C>7 
is probably determined by a delicate balance between 
the long-range dipolar interaction and several weak next- 
neighbor exchanges^ The full analysis of the role of the 
single-ion anisotropy for such models lies beyond the 
scope of present study. We will discuss instead properties 
of a simplified classical spin model on a pyrochlore lat- 
tice, which takes into account only the nearest-neighbor 
exchange interaction and a staggered single-ion term: 

H = J^Si-Sj + D^fa-Sif-H-^Si ■ (3) 

(ij) » i 

Transformation to classical spins of unit length |Sj| = 1 
amounts to the following redefinition of the physical pa- 
rameters: J -> JS 2 , D DS 2 , H -> g(J, B HS. The 
minimum of the Heisenberg interaction in zero field is 
attained for states with zero total spin on every tetrahe- 
dron: S tc t = Yli=i &i = 0i£ The key observation for the 
anisotropic model J3J) is that the Stct = constraint is 
compatible with local spins being oriented perpendicular 
to their local hard axes: the exchange and the anisotropy 
terms can be simultaneously minimized. 15 -'-- In partic- 
ular, for one tetrahedron there is one free continuous 
parameter (angle), which describes the possible ground 
states. The remaining degeneracy of the lattice model 
is infinite but not extensive: the number of continuous 
parameters describing the classical ground state scales 
as 0(L 2 ) with the linear system size L. Monte Carlo 
investigations^ 5 ^ of the model © with D/J — 5 and 
00 have found a first order transition into a noncopla- 
nar L q — 0' state at T c ~ 0.13J, which is driven by the 
thermal order by disorder mechanism. 

To complement the above zero-field consideration we 
now discuss the magnetization process of the model (J3J- 
The total spin operator does not commute with the 
Hamiltonian in the presence of anisotropic interactions. 
Consequently, the magnetization plateaus, e.g., the 1/2- 
plateau, do not appear on the magnetization curves of the 
gadolinium pyrochlores and M(H) approaches the satu- 
ration value at high fields only asymptotically. Besides, 
the staggered (within one unit cell) single ion anisotropy 
might also smear phase transitions in finite magnetic 
fields for 'nonsymmetric' field orientations. 

Returning back to the Hamiltonian the minimum 
of the Heisenberg and the Zeeman interactions is given 
by the states obeying Stct = H/2 JJ^ Simple extension of 
the zero-field consideration suggests that in weak mag- 
netic fields spins remain in their respective local easy- 
planes. It is again possible to simultaneously minimize 
the local single-ion anisotropy terms and to satisfy the 
magnetization constraint. The magnetization curve is 
initially a straight line: M{H) = H/8J. The remaining 
degeneracy corresponds to one continuous degree of free- 
dom per tetrahedron and 0(L 2 ) variables for the lattice 
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FIG. 3: (Color online) Magnetization per spin for the clas- 
sical easy-plane pyrochlore antiferromagnet at T = for the 
field applied along the two symmetry directions. Magnetiza- 
tion curves from top to bottom are drawn for D/J = 0.2, 1, 2, 
and 5. The transverse magnetization on the right panel cor- 
responds to the same choice of D/J in the reverse order. 

problem. 

The situation, however, changes with increasing mag- 
netic field. At a certain critical field H c , which is 
smaller than the saturation field of the Heisenberg model 
H sa ,t = 8JS, the constraint Stct = H/2 J cannot be sat- 
isfied anymore by spins lying in their local easy-planes. 
The degeneracy is lost and spins start to tilt out of the 
easy-planes towards the field direction. The correspond- 
ing critical field does not depend on the value of the 
anisotropy constant and is determined purely by geom- 
etry It can be calculated analytically for fields applied 
along the symmetry directions of a pyrochlore lattice: 

tfi 0011 = 8 J ^2/3 , ff| lu l = 8J(3 + V2)/7 , 
ffl 1121 = 6 J , Hi 1101 = 4J(1 + l/\/3) . (4) 

The fields 7?I 001 ' and if! 1 give, respectively, the maxi- 
mal and the minimal value for H c (8, <p). The magnetiza- 
tion process at T = for the above two orientations found 
by standard numerical minimization of the classical en- 
ergy of a single spin tetrahedron is presented in Fig. 
For H || [001] the magnetization curve exhibits a kink 
at H c and, then, the behavior of M(H) becomes non- 
linear with full saturation M = 1 reached only asymp- 
totically. A nondegenerate spin structure above H c pre- 
serves the C2 rotational symmetry about the field direc- 
tion, which is broken at H < H c . Thus, H c describes 
a robust phase transition for this field orientation. The 
zero-temperature magnetization curve for arbitrary field 
orientation resembles the behavior for H || [001]. The 
main difference is an absence of a spin-rotational symme- 
try about the field direction, which may lead to a smear- 
ing of the phase transition at H — H c . The other two 
field orientations, which have stable phase transitions at 
H = H c are those for the field along the [110] and the 
[111] axes. In the latter case the magnetization curve ex- 
hibits two kinks: at H c \ = H c k, 5.045 J (0} and at a sec- 
ond critical field H C 2 > H sat , see the right panel in Fig. 
Above H C 2 the spin structure preserves the C% v symme- 
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FIG. 4: (Color online) Transverse magnetization per one spin 
of the classical easy-plane pyrochlore antiferromagnet at T = 
for D = 0.75J and various field orientations. 

try of the system in magnetic field, while for H < H C 2 
only a mirror symmetry remains with respect to one of 
the three mirror planes passing through the [111] axis. 
Such a remaining symmetry is further lost in the low- 
field phase at H < H c \. The upper transition field has 
an almost linear dependence on the anisotropy constant 
in a wide range of values for D < 10: H C 2 ~ 8 J + 1.7D. 
We, therefore, conclude, that for H || [111] there are two 
phase transitions in magnetic field, which are smeared 
when the field is tilted away from this symmetry direc- 
tion. 

In the intermediate phase H c \ < H < H C 2 with 
broken three- fold rotational symmetry (H || [111]) the 
anisotropic pyrochlore antiferromagnet develops a trans- 
verse magnetization, which is shown in Fig. [3] The ampli- 
tude of the transverse magnetization grows quickly with 
increasing anisotropy constant D. There arc, of course, 
three domains of Mj_ oriented under 120° to each other 
in accordance with a broken symmetry of intermediate 
phase. Once magnetic field is applied away from the [111] 



axis, Mi attains a unique direction. To make connection 
with Gd2Ti2C>7, which according to our measurements 
has D « 0.75 J, we plot in Fig. 0] the amplitudes of the 
transverse magnetization for the above ratio of D / J and 
for several characteristic orientations of the applied field. 
Though Mj_ does not exceed 3% of the saturated magne- 
tization, the corresponding values can be definitely mea- 
sured in an experiment. Such measurements would also 
serve as a precise method to determine the two phase 
transitions for H j| [111]. The transverse magnetiza- 
tion has been previously measured in CsMnBr3, a tri- 
angular lattice antiferromagnet with a strong single-ion 
anisotropypi& Note, that for H || [001] and [110] the trans- 
verse magnetization does not appear in the whole range 
of magnetic fields. 

In conclusion, the measured local anisotropy of Gd 3+ 
ions strongly affects magnetic properties of Gd-bascd 
pyrochlorcs. Though, the toy model J2J) is apparently 
not sufficient to explain the ordered structure seen by 
neutrons^ the quite anisotropic phase diagram in mag- 
netic fields is compatible with our experimental and the- 
oretical findings. In particular, we have demonstrated 
that depending on an applied field orientation in a cubic 
crystal there are possible cither one or two transitions 
in magnetic field. The Mossbauer measurements^ also 
agree with our ESR results pointing at ordered moments 
oriented perpendicular to local trigonal axes. The pre- 
sented theory of the magnetization process is relevant for 
Er2Ti207, which is a D/J — > oo realization of the model 
©■ 
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